Contribution of Light-by-Light Scatteringto Energy Levels of Light Muonic Atoms 
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The complete contribution of diagrams with the hght-by-hght scattering to the Lamb shift is 
found for muonic hydrogen, deuterium and helium ion. The results are obtained in the static muon 
approximation and a part of the paper is devoted to the verification of this approximation and 
analysis of its uncertainty. 
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Studies of energy levels in muonic hydrogen have a long 
history, but until recently there have been no successful 
precision measurements on this atom. Recently results 
on the Lamb shift in muonic hydrogen and deuterium 
were obtained for the first time [l| , and similar measure- 
ments for muonic helium are planned. 

The Lamb shift in light muonic two-body atoms is a 
splitting of the levels 2s and 2p, which are degenerate 
within a nonrelativistic (NR) treatment of the Coulomb 
problem (as well as, e.g., in common hydrogen). Ef- 
fects of the electron vacuum polarization in muonic atoms 
break down the degeneration and lead to the splitting of 
order a{Za)^m^, which can still be obtained from NR 
calculations. (The relativistic units, in which h = c = 1, 
Z is the nuclear charge and a is the fine structure con- 
stant, are applied throughout the paper.) 

There are also "finer" splittings and, in particular, lev- 
els in light muonic atoms possess the fine and hyperfine 
structure, which result from relativistic calculations and 
appear to be substantially smaller than the Lamb shift 
(in contrast to common hydrogen). 

Specificity of light muonic atoms is that the charac- 
teristic atomic momentum is Zam^, comparable to the 
mass of electron (am^ ~ I.SttIe). At the same time, 
atomic energies are much lower than me- 

Therefore, one can consider corrections to energy with 
closed electron loops as nonrelativistic. Effects of the 
vacuum polarization up to NR contributions to the Lamb 
shift of order a^m^ were studied in a number of papers 
(see, e.g., (H). 

While calculating a^m^ terms in the papers men- 
tioned, the atomic nucleus was treated as a point-like 
one. Nuclear-finite-size effects in light muonic atoms can 
be considered, if necessary, as an additional perturbation. 

Contributions, induced by the light-by-light (LbL) 
scattering (see Fig. [T]), appear in order of a^m^ and are 
considered in this paper. The atomic nucleus in diagrams 
in Fig. [1] is treated in the external field approximation 
(i.e., as a static one), while the expression for the muonic 
line includes the Coulomb Green function of muon, which 
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FIG. 1: Three types of diagrams for the LbL contributions to 
the Lamb shift in light muonic atoms 



is indicated by the double line. 
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The sum is taken over all intermediate states A of a dis- 
crete and continuous spectrum. 

Each of the photons, including those that connect the 
muon line and the electron loop as well as those of the 
external nuclear field, are Coulomb photons, i.e., only 
the Z3qo component of the photon propagator contributes 
to the result. Z?oo does not depend on energy in the 
Coulomb gauge, however, energy is transferred through 
photons lines that connect the muon line and the electron 
loop. Apparently, energy does not propagate through the 
external field lines. 

Generally speaking, since the muon atomic momentum 
is of order of the electron loop momenta, explicit forms 
of the functions |A(p)) are substantially different for free 
and bound (Coulomb) wave functions (see, e.g., 0). 

Calculations with Coulomb Green functions turn out 
to be rather complicated even in the NR approximation; 
up to date they have not been performed for any of the 
corresponding contributions. 

However, it is possible to demonstrate that the calcula- 
tions can be made in a simple approximation considering 
muon as a static particle (see Fig. [5]). 

In this approximation the expressions for individual 
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FIG. 2: Diagrams for the calculation of the Lamb shift con- 
tribution in the static muon approximation 



contributions to the energy of the level C, are of the form 
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where 



^d^^=J^nip)^dp 



(5) 



is the atomic form factor, which is for the 2p—2s splitting 
found to be 



J^2p-2s{ci) 
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7 = ZaTTir, and tti^ is the reduced mass of the muon in 
the atom. 

The factor £(fci, A;2, fca, ^4) corresponds to a single di- 
agram of the LbL scattering with the momenta of the 
incoming photons defined as ki = (0,ki) for i = 1,2,3 
and k4 = — (fci + A;2 + fca) (see Fig. [3]). The permuta- 
tions of photons, as it is demonstrated below, have been 
already taken into account explicitly. 

We remark that both the Coulomb Green functions 
and integrations over photon energy that are included 
in the full expressions have completely vanished. Note 
that energy in the static approximation does not prop- 
agate through photon lines, but does through photons 
that connect the muon and electron. 

The static muon approximation qualitatively changes 
the form of expressions and notably simplifies the calcu- 
lations. 



FIG. 3: The diagram for the LbL scattering block 
C{ki,k2,k3,k4) 



We present here the proof of its applicability to light 
muonic atoms, assuming at the first stage that the bound 
energy of the state C, for which the energy is calculated, 
as well as for intermediate states A, is negligible. Veri- 
fication of this assumption is discussed separately after- 
wards. 

The standard muonic Coulomb Green function appears 
in our calculations with the following arguments; 

SciEc±u:,p,p'), 

where w is a photon frequency or combination of such 
values. In contrast to the static muon approximation 
(see Fig. [2|), these frequencies are not zeros in the full 
expression. 

Neglecting the energies of the states and taking into 
account completeness of the basis, one can find that the 
sum over A becomes trivial: 
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Expressions with the Coulomb Green function differ 
from expressions in the static muon approximation be- 
cause three additional elements appear: 

1) a sum over intermediate states; 

2) an integration over momentum (note, while a muon 
propagates at the Coulomb field, its energy rather than 
momentum is conserved); 

3) an integration over the photon frequencies (note that 
the virtual photons, connecting the muon and electron 
in Fig. [1] transfer energy, while the photons in the static 
approximation in Fig. [3] do not). 

One should also take into account contributions of all 
possible permutations of the photon lines. We treat in 
our calculations a single diagram in Fig. [3] as the whole 
LbL scattering block, and all the permutations are taken 
into account as contributions to the muonic factor with 
permuted lines of outgoing photons. 

The transformation (JT)) allows removing both the sum 
over intermediate states and "superfiuous" integrations 
over momentum. Now the only difference between the 
expressions ([3]) and ^ and the Coulomb one is the in- 
tegration over energy of the photons. (The expression 
(12) for the l:3-contribution is already equal to the full 
expression since it does not really include the Coulomb 
Green function from the beginning.) 
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FIG. 4: The muon factor for the 2:2 contribution 



Let us look at the structure of integrations over u. 
Exchange of energy between the muon and electron is 
present in the Coulomb diagrams, but complete energy 
transfer from the nucleus to electron is equal to zero, and 
therefore the sum of energies of all photons connecting 
the electron and muon is zero as well. Thus, in the case 
of the 2:2 contribution, having defined the energy of the 
first photon fcio as oj (so that fc2o = — w), we find that 
the sum of the direct and cross diagrams for the muon 
factor (see Fig. ^ leads to a 5 function 



oj + iO —u! + iO 



= 27r(5(w) 



(8) 



that removes the last "superfluous" integration. As a 
result, we obtain the expression (|3]). 
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FIG. 5: The muon factor for the 3:1 contribution 

The structure of the 3:1 diagrams is more complicated. 
Having deflned energies of the photons as kio — wi, 
feo = ^2 (so that fcso = — (wi + W2)), the sum of six 
permutations for the total muon factor (see Fig. [5]) takes 
the form 
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The 5 functions remove integrations over energy and lead 
to the static expression Q. 

Therefore, neglecting the energy of the states in the 
Coulomb Green function ([71), we have in fact turned the 
Coulomb diagrams in Fig. [T] into diagrams in Fig. [5] with 
a static muon. We still have to verify that it is really 
possible to neglect the dismissed terms and estimate the 
error of this procedure. 

In case of the 2:2-contribution the structure of the in- 
tegral over the photon energy for individual diagrams is 
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where we define the combination of energies that we are 
to neglect as e ^ mg. 

The first factor is for the muon and it is distinct for 
different contributions and the second one, C{uj)^ which 
corresponds to the LbL scattering block, is universal. 

After combining with other diagrams the first factor 
leads to the 5 function at zero frequency (cf. © and 
dH])), therefore, strictly speaking, it cannot be expanded. 
One can see that after performing a substitution of the 
variable by shifting w — )> a;' =p e, the integral turns to 
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This expression can already be expanded in powers of e. 
It is clear that while integrating in the LbL block, the 
characteristic scale of loop energies in C{uj) is of order 
of the electron mass (or higher) and the parameter of 
expansion is 

±_ ^ ^Zafm, ^ ^^^ ^2 
nif, nif. 

It is clear that a similar approach with substitution 
of a variable is applicable for the 3:1 correction as well. 
E.g., for the first term in ([9]) one can write 
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and then the expansion in e can be done without prob- 
lems. 

We note that, if the diagrams in Fig.[T]contain the free 
muon propagators instead of the Coulomb Green func- 
tions, the muon static approximation is applicable and 
the uncertainty is of the same order. Dealing with the 
the muon Green function (see, e.g., ([7])), only the com- 
pleteness of the eigenstate basis is required, while the 
estimated eigenvalues of the energy are of the same or- 
der for the free and Coulomb case, as long as the related 
integrals are convergent. 

One of direct consequences of the applicability of the 
static muon approximation and of s^jonmetry of the ex- 
pressions ^ and (JH) is the identity 



Ai?3:i — -—^AEi;3 



(12) 



Until recently the 3:1 term has been the only com- 
pletely unknown contribution. On the contrary, the 1:3- 
contribution, which is also referred to as the Wichmann- 
KroU contribution, has been well known (see @-|g| for 
results in muonic hydrogen, [9| in muonic deuterium and 



[3, [lO] muonic helium). That is due to the fact that 
approximations of the Wichmann-KroU potential in the 
form of an explicit function of distance, which are effi- 
cient for calculations, have been known [ll|, [ij]- There 
is also an exact representation of this potential in the 
form of a double integral [l^] , understood in terms of the 
principal value. 

The identity p^ changes the situation radically. It 
by itself improves accuracy of contributions to the Lamb 
shift in muonic hydrogen and muonic deuterium by an 
order of magnitude and also leads to a certain improve- 
ment of the accuracy for muonic helium. 

According to Eq. p^ . the 2:2 term, which is also re- 
ferred to as the virtual Delbriick scattering, becomes the 
least accurately known contribution from the LbL block. 
The results for muonic hydrogen |14| and deuterium (9[ 
are obtained with accuracy less than 10%, and for muonic 
helium-4 ion it has been rather an estimate [10| than a 
result obtained. 

The 2:2 contribution was examined in the cited papers 
in the so called scattering approximation, where at first 
an operator related to the Feynman diagram with exter- 
nal muon lines on the mass shell and with a free muon 
propagator, is derived, and then its matrix element over 
the Coulomb wave function is calculated. In fact, before 
the calculation was really done [ij, \1^ , a few additional 
approximations had been introduced, which reduced the 
applied expressions to the static muon approximation. 

Leaving aside the question about reasonability of using 
the scattering approximation as an initial point, we state 
that direct calculations in [3, Hfl Il3| were in fact per- 
formed in the static muon approximation, which is cor- 
rect for the corresponding Coulomb diagrams in Fig. [T] 
within the declared accuracy of calculations. 

The 2:2 contribution is one of the most complex spe- 
cific QED contributions in light muonic atoms and ab- 
sence of independent confirmations significantly reduces 
reliability of the results obtained in [3, \l^ [iJI . 

Below we calculate the 1:3 and 2:2 terms. The first 
calculation is used for control of analytical expressions 
and numerical algorithms. The second one is aimed to 
improve accuracy of the total contribution of the diagram 
in Fig. [T] for light muonic atoms, since the accuracy of 
known results is not high enough. 

To perform a calculation, one has at first to find an effi- 
cient form of the factor C, which corresponds to the LbL 
scattering diagram (see Fig.[3|). A calculation within the 
standard technique of Feynman parameterization leads 
to the result 
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Here, we have explicitly introduced the Pauli-Villars 
regularization with M ^ Trie, which can be useful if the 
convergence of the integrals over k is not good enough. 

It is useful to make a further calculation by combining 
the photon denominators and the denominator A using 
Feynman parameters. After integration over momenta 
we come to the following expression: 
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and the dimensionless coefficients Ai.s, Bi;3, Ci:3 and si:3 
are bulky functions of all Feynman parameters. 

Calculating the 1:3 contribution, the integration over 
the momentum of the atomic form factor J^ in ([2]) is fac- 
torized, and the remaining integrations over momenta in- 
volve logarithmic divergencies at large momenta. There- 
fore, keeping M at intermediate stages turns to be useful 
for a calculation of individual terms. In particular, con- 
tributions to (|14p of the separate terms of P^ include 
such divergences. 

The value J^i-.siq^) is nothing else but a contribution to 
the charge form factor of the muon induced by the LbL 



block. We can always "renormalize" the vertex function 
of muon and, by subtracting 



Jl:3(q')^ Jl:3(q')-Jl:3(0), 



(16) 



remove the logarithmic divergency in integration over k 
in (HSl). 

In fact, the LbL scattering diagram, as it is known, 
does not renormalize the vertex, i.e., J7i:3(0) = 0, that 
we have checked both numerically and analytically at 
different stages of transformations. In particular, the 
coefficient at the logarithm ln(M/me) in ([T5|) turns to 



dy dz I du I dv dw Ai-. 



dx 



Therefore "renormalization" of individual terms accord- 
ing to P^ is actually reduced to regrouping and cance- 
lation of individual divergent terms. 

The final integration over Feynman parameters and 
momentum q has been performed by numerical integra- 
tion with VEGAS Ta] . The results are presented in Ta- 
ble m The contributions for the 2p — 2s splitting in the 
last line of the table were calculated directly rather than 
as the difference of the 2p and 2s terms. 



Level(s) C(/iH) C(mD) 



CifiUe) 



2s 0.5704(5) 

2p 0.10468(3) 

2p-2s -0.4649(5) 



0.6263(5) 1.0815(10) 

0.12405(3) 0.5108(7) 

-0.5015(5) -0.5702(10) 



TABLE I: The 1:3 contribution (the Wichmann-Kroll contri- 
bution) to energy levels of muonic hydrogen, deuterium and 
helium-4 ion: AE = a{Za)'^mf, ■ 10"^ ■ C 



The results for the Lamb shift 2p — 2s (see also 
Table IIIip are in excellent agreement with the val- 
ues previously obtained by other authors for muonic 
hydrogen (-0.00103(2) meV [I, i]) and deuterium 
(—0.00111(2) meV 0). For helium ion our result agrees 
with -0.0198(4) meV 0] and -0.02 meV [H, but dis- 
agrees with the result 0.135 meV [l^]- The uncertainty 
of the results [3, Q was not present in the original papers 
and was estimated here according to Q. 

The momentum integrations for the 2:2 contribution 
were performed similarly to those described above. One 
can check that it is possible here to set M = oo be- 
fore integration over large momenta. After introducing 
Feynman parameters and integrating over momenta, we 
obtain 
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•^2-2' Q-2 ' '^2-2 ^'^d 4-2 ^'"^ bulky dimen- 
sionless functions of the Feynman parameters. Note that 
the denominators for the ladder (fc = 1) and non-ladder 
{k = 2) diagrams are slightly different. 

The integration over q and the Feynman parameters 
is performed numerically. The results are presented in 
Table HIl Similarly to the case of the 1:3 contribution, 
the value for the 2p — 2s splitting in the table is found 
by a direct numerical integration. 



Level(s) CjfiR) 



g(MD) 



CifiUe) 



2s -0.8201(15) 
2p -0.2937(9) 
2p - 2s 0.5264(12) 



-0.9000(15) -1.5615(25) 
-0.3355(9) -0.908(2) 
0.5645(12) 0.652(2) 



TABLE II: The 2:2 contribution (the virtual Delbriick scat- 
tering) to energy levels of muonic hydrogen, deuterium and 
heIium-4 ion in the static muon approximation: AE = 



a^iZafrnf, • lO""^ -C 



The results for 2p — 2s in the static muon ap- 
proximation are slightly smaller and have substantially 
higher accuracy than the results of the other authors. 
Ours are mostly in agreement with the former val- 
ues. In particular, for muonic hydrogen we obtain 
0.001 151(4) meV (cf. 0.001 35(15) meV (H), for deu- 
terium, 0.001234(4) meV (cf. 0.00147(16) meVjl) and 
for hchuni ion, 0.01140(4) meV (cf. 0.02 meV fiof). 

Complete contributions of the diagram in Fig.[T]to the 
Lamb shift are collected in Table Hill In contrast to Ta- 
ble nil where the uncertainty of numerical calculations in 
the static muon approximation is presented, here we take 
into account the error of the static muon approximation, 
which dominates. 

In conclusion, we note smallness of the numerical co- 
efficients in Tables U and HIl which is typical for spin- 
independent nonrelativistic contributions of the light-by- 
light scattering to energy levels. 

As we have noted, the results for the 2:2 contribution 
obtained in this paper are slightly below the results of 
Boric [3, \1m> uM- The difference for muonic hydrogen 
and deuterium is about 20%, i.e. less than 1.5 stan- 
dard deviations. In the case of muonic helium our re- 
sult (0.01140(2) meV) is also lower than that of Borie 
(0.02 meV), however, the latter is presented in [10| in 



Term 



AE{fj.R) 

[meV] 



AEinB) 

[meV] 



A£(/xHe) 

[meV] 



1:3 -0.001018(4) 
2:2 0.00115(1) 
3:1 -0.00102(1) 



-0.001 098(4) 

0.00124(1) 

-0.00110(1) 



-0.019 95(6) 

0.0114(4) 

-0.0050(2) 



Total -0.000 89(2) -0.000 96(2) -0.0136(6) 

TABLE III: Contribution of the LbL scattering effects 
(Fig. [l]) to the Lamb shift (2p — 2s) in muonic hydrogen, 
deuterium and heUum-4 ion 



such a form that one is rather able to estimate the scale 
of uncertainty than its magnitude. 

This difference in the results can in no way be consid- 
ered as a contradiction, but it has a systematic character 
(all our results are lower than those of Borie) and, strictly 
speaking, requires some explanation. 

Due to that, we have checked two most important ele- 
ments of our calculations, namely, the expression for the 



LbL scattering (TTSl) and the integration over momenta 
ki in (U]). As a test of the representation for £ we have 
used a calculation of the well-known 1:3 contribution (see 
Table IJ) and studied the expression (IT^ for the LbL scat- 
tering in the case when momentum of one of the photons 
is zero. 

To be certain that there are no systematic errors in 
integration over momenta after introducing Feynman pa- 
rameters we have made an independent calculation of the 
2:2 contribution, in which integrations over momenta ki 
in the corresponding expression (J4l) were performed di- 
rectly. The results turn out to have accuracy 4% for 
muonic hydrogen and 3% for muonic helium ion and are 
in excellent agreement with results of Table |lll 

These tests allow to consider our results to be highly 
reliable. 
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